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Abstract 

We show that a spontaneous bending of single layer bolaamphiphiles results from the frustration due to 
the competition between core-core and tail-solvent interactions. We find that spherical vesicles are stable 
under rather general assumptions on these interactions described within the Flory-Huggins theory. We con- 
sider also the deformation of the vesicles in an external magnetic field that has been recently experimentally 
observed. 

Molecular aggregates, like biological matter, possess the ability to organize into well defined mesoscopic 
structures, such as layers, fibers, vesicles PEl El E]- Addressing this issue for molecular aggregates is relevant not 
only for applications in pharmacology, catalysis, and other fields, but also to understand the behaviour of model 
systems simpler but related to biological matter. The self-organized structures result from the concomitant 
effect of non-covalent interactions, such as hydrogen bonding and ir — ir stacking between molecules, as well as 
by entropic contributions. A theory capable of predicting the shape of the aggregate from the knowledge of the 
microscopic molecular structure is still missing. Among the forms of self-organization, the formation of empty 
vesicles in, possibly aqueous, solution is particularly important for the wealth of possible applications, from 
microreactors to drug delivery. Surprisingly, vesicle formation has been observed not only in the well known 
case of amphiphiles bilayers, like lipid membranes but also in several bolaamphiphiles[Hl[I|IH] formed by 
an hydrophobic core and two, usually symmetric, lateral hydrocarbon tails terminated by hydrophilic groups|Hj- 
Notice that the term bolaamphiphiles has been first introduced to describe synthetic analogs of archaebacterial 
membranes |lUI ITT] . Here we provide a model predicting spontaneous bending due to frustration resulting from 
competing core-core and tail-solvent interactions. 

We refer, in particular, to the case of sexithiophenes bolaamphiphiles (see Fig. that have been found to 
form vesicles albeit in isopropanoljHj- The formation of stable vesicles seems to be favoured by increasing length 
of the tails jjj. Conversely, layer- like structures have been found in liquid crystal bolaamphiphiles with short or 
rigid lateral blocks [T2l 113) . 





Figure 1: Structural formula of 2,5"" '-(R- 2- methyl- 3,6,9, 12, 15-pentaoxahexadecyl ester) sexithiophene and its 
schematic representation. 

We consider a simplified model of a bolaamphiphile layer with the cores on a two-dimensional surface and 
symmetric flexible tails, each formed by m monomers as sketched in FigEI 

As illustrated in FigEl the cores define a surface So with local coordinates (u\, U2) such that r*o = r(ui, U2) 
gives cartesian coordinates on the surface. If the tails, due to covalent bonding to the cores, tend to be 
perpendicular to the surface, the interactions with the solvent are supposed to take place on the two surfaces 
S± with coordinates 

r± = f(ui,u 2 ) ±n(ui,u 2 )D (1) 

where n is the vector normal to So and D is the effective length of the tails. 
The elementary area vector on So is 

dSo = {r\ x r%) du\du2 (2) 
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Figure 2: Sketch of a layer formed by self-assembled bolaamphiphiles. The spacing between hydrophilic tails (green) is 
fixed by the hydrophobic cores (yellow) to a value which might be not optimal for tail-solvent interactions. 



where = dr/dui. The derivative of the normal n with respect to Ui is given by the so-called Weingarten 
equations [21, whence 

dS± = dS (l =F 2HD + KB 2 ) (3) 

where H = («i + K2)/2 and K = k±K2 are, respectively, the mean and gaussian curvature defined in terms 
of the principal curvatures Ki and K2- The tail density in contact with the solvent is determined by the tail 
density no imposed by the core so that n±dS± — ngdSo leading to an exact expression for the tail density on 
S± in terms of the curvature 

710 (A\ 

n± = 4) 

1 =p 2H D + KB 2 v ; 

If the tails were not rigidly connected to the cores, they would form a solution with free energy per unit volume 



Figure 3: Schematic geometry of our model. The surface So is the locus of the center of mass of the hydrophobic cores 
and S± are the surfaces formed by the tail ends. The unit vector n is normal to all three surfaces and D is the effective 
tail length. 




4>(n). Therefore, we write the bending free energy as 

5F = D J dS [<t>{n+) + 0(n_)] + A J dSH 2 +B J dSK. (5) 

The first term is the free energy of the tail/solvent subsystem whereas the terms proportional to A and B 
represent the bending elasticity of the core layers El- To account phenomenologically for bent surfaces 
for lipid bilayers, the term AH 2 is usually replaced by the so-called Canham-Helfrich (CH) bending energy 
A(H — Ho) 2 with A > so that negative linear terms in H are responsible for bending^. However, the CH 
expression is not appropriate for symmetric bolaamphiphiles where odd terms vanish by symmetry. Moreover, 
the spontaneous curvature Hq is introduced phenomenologically and is not justified by any microscopic model. 
We also note that the term B J dSK is usually neglected because due to the Gauss-Bonnet theorem^ a it is 
constant within a given topology of the surface and it actually depends on the number of handles. We consider 
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mesoscopic vesicles with curvature much larger than the molecular size, namely \Ki\D << 1. Within this 
approximation we can expand <b(n±) up to second order in D as 

<t>(n±) = 0(no) ± n Q( j)'(n )2HD+ 



+no4>'(no) [4H 2 D 2 - KD 2 } + ^nfo" \n )AH 2 D 2 (6) 

where <fi'(no) — d(bjdn\ n=na and <f>"(no) = d 2 (b/dn 2 \ n=no . Notice that the sum of the terms linear in H in the 
integrand of Eq[S] vanishes whereas other terms contribute to the term proportional to H 2 so that 

A—>A = A + 5A = A + 4£> 3 [nl<t>" + 2n </>'] (7) 

We will show that A can become negative and lead to spontaneous bending as a result of the microscopic 
tail/solvent interactions. The reason for a negative 5A is that the available volume per tail on the two sides of 
a bent surface may be different and a specific curvature can minimize the resulting free energy. 

We can further specify the driving term 5 A by calculating it, within the Flory-Huggins theory for polymers 
solvent mixtures |14l IT] 

(bin) mnni 
1 rii + mn 

+ n 1 ln(^-)+n 2 ln(^^-) (8) 



rii + mn J \ni 

where n\ is the solvent density, to is the number of monomers per tail, \ ~ (2Wt« — W ss — Wu) /T with Wt s , W ss 
and Wu the tail-solvent, solvent-solvent and tail-tail interaction energies respectively and T is the temperature. 
Substitution of Eq.8 into Eq.7 gives 



T (m + mno) 2 (m + mno^ 2 



• [m (to — 2) no + (2to — 3) m] (9) 

For hydrophilic tails, x < leading to a free energy gain for bent surfaces. This constitutes a possible micro- 
scopic mechanism for vesicle formation, albeit within the Flory-Huggins theory which is more qualitative than 
quantitative [T^]. Furthermore, the last term that accurately represents entropy is always negative unless m = 1 
and increases as to 2 for to 3> 1. The observation of planar structures in liquid crystals bolaamphiphiles with 
rigid tails could be explained by this entropic reason. This completes our microscopic analysis. Our model 
leads naturally to a negative contribution to the free energy of the term quadratic in the curvature that does 
not need to be imposed phenomenologically as done for amphiphile bilayers. 

Our results can be cast in the form of an effective Landau Hamiltonian, i.e. we can write the free energy F 
as a Taylor expansion in terms of the curvatures H and K. These quantities, or more precisely K\ and take 
the apparent role of order parameters assuming their smallness in comparison to a typical inverse molecular 
size. However, n\ and n 2 are not independent at different points of the surface, implying that F is not a true 
Landau Hamiltonian and that a local analysis is insufficient. In fact, the local free energy density has not a 
minimum but a saddle point for a sphere (k% — k 2 ) if c 2 + 2c% < 0, in contradiction with the numerical data 
presented below. Since the term proportional to H 2 resulting from our microscopic model is negative, higher 
order terms are needed to stabilize the mean curvature to a finite value. Up to fourth order terms in k\ and k 2 
we can write 

F = J dS[-\A\H 2 + bHK + b'H 3 + c x H A + c 2 H 2 K + c 3 K 2 ] (10) 

The odd terms in b and b' that can lead to first order phase transitions have to vanish for the symmetric 
bolaamphiphiles considered here so that we take 6 = and b' = 0. Moreover, c\ + c 2 + C3 > 0. Since, in general, 
the bending energy is much smaller that the surface tension, we go to the limit of infinite surface tension that 
allows to perform minimization of F on a surface of constant area. We have minimized Ea llOl numerically for 
surfaces topologically equivalent to spheres with radius R = y/2 (c\ + c 2 + C3) /\A\ and found that spheres are 
always the solution. 
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Our model can be extended to consider the recently observed effect of external magnetic fields B on the 
equilibrium shape of bolaamphiphilic vesicles To this purpose we add to Eg 1101 the diamagnetic energy term 

E mag = ~ £ J dS ( X i sin 2 9 + X \ cos 2 6) (11) 

where Di is the effective thickness of the core layer (t = 1) or of the tails (i = 2), Xj_j X\ are the components 
of the corresponding magnetic susceptibility and 9 is the polar angle of the direction of the magnetic field 
B. Due to the tendency of the aromatic rings in the core segment to align parallel to the magnetic field, 
the spherical vesicles will be transformed into ellipsoids, as predicted by Helfrich|17| for amphiphilic bilayers. 
This deformation in high magnetic fields has been recently observed experimentally for sexithiophene 8 , that 
is an example of symmetric bolaamphiphilcs 8 where the CH model of free energy is not correct. In FigQ] 
we present the deformation calculated by minimizing the sum of Eqs. 10 and 11 for constant surface as a 
function of magnetic field, compared to the one resulting from the CH model. In both cases, the deformation is 
proportional to B 2 for low fields and flattens out as the field increases. However, for the same initial slope, our 
model predicts smaller deformations at high fields, the details of the curve being determined by the parameters 
ci, c 2 , c 3 . 
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Figure 4: Difference AR between short and long axis of the ellipsoid resulting from minimization of the free energy 
in a magnetic field B for three choices of C2 and C3 (ci + C2 + C3 = 1 and \A\ — 1), compared with the Heffrich model 
(Ho = 1). The inset shows that initial B 2 dependence of the distortion depends mostly on C2 + 2c A . 

In summary, we have demonstrated that the frustration resulting from competing core-core and tail-solvent 
interactions can lead to spontaneous bending of single bolaamphiphilic layers. By describing the tail solvent 
interactions within the Flory-Huggins theory, we have constructed a Landau-like free energy appropriate to 
describe symmetric bolaamphiphilcs that gives a rationale for the formation of spherical vesicles. Measurement 
of their deformation in high magnetic fields can provide information about the parameters of the theory, opening 
the possibility to validate microscopic models of interactions in these systems. 
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